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The classical walking behaviors of a single atom in an amplitude-modulated standing wave lattice 
beyond the internal dynamics are investigated. Based on a simple effective model, we identify a 
diversity of dynamic regimes of atomic motion by periodically adjusting the lattice depth. Harmonic 
oscillation or pendulum rotation with classical step-jumping, random scattering walking, chaotic 
transportation, quasi-periodic trapped motion and roughly ballistic free flying are found in this 
simple model within different parametric regions by approximate analyses. Our study demonstrates 
' a complex motion of single atom in modulating optical lattice beyond the quantum description. 
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I. INTRODUCTION 



Nowadays, the fabrication technology steps more and more towards the atomic scale, and detection on one single 
molecular or atom is now an promising technique 0, [|| . Controlling one atom in the cavity quantum electrodynamics 
has now nearly met this step [3H5J . Based on this situation, the motion of one single atom in a controlled 1-dimensional 
(ID) optical lattice becomes a hot topic with a continuous interestingfor its basic theoretical value and reliable 
applications in quantum information processing Q , atom lithography [7| , solid state physics Q , quantum walking 
algorithm || [l(|, plasmonic application (ill [Tpl . quantum-to-classical transition [l3| and molecular dynamics. In the 
optical lattice, the atom can be accelerated |l4j or decelerated [l5[ by an optical force for the purpose of precise 
control. However, in atomic scale, quantum effect starts to play a tricky role in the atomic dynamics. For example, 
under a deterministic linear force, the atomic motion is still beyond definite prediction in quantum theory because of 
^ ' an intrinsic uncertain variation. That the quantum motion distinguished from a classical trajectory motion is actually 
due to the superposition of the atomic wave packet which induces interference effect such as dynamical localization 
But, just for one single atom, the external wave superposition will be excluded [l7j and the motion at different 
time, also, will not interfere unless a long memory imprinted on the optical field can feed back to the atomic motion 
along the classical trajectory. However this memory will also be cut down by a definite external optical control. 
Therefore the quantum motion and the classical motion meet together at this single atomic scale. 

According to quantum descriptions by including internal level transitions, many kinds of motional behaviors were 
found in this single atomic system, such as nonclassical motional states, ballistic transport, harmonic oscillations, 
random walks, Levy flights, and chaotic transport [1 81 - 4231 ] - Which of these behaviors are derived from quantum 
description and which comes from classical dynamics are still a mixing problem in this system. This multiple dynamics 
inspires us to find if there exists a simple classical model without internal dynamics to reproduce most of the above 
dynamical processes. The interesting thing is that, in this paper, we do find a simple classical model which can 
produce most of above behaviors without directly considering internal dynamics. We identify a rich dynamic picture 
of the atomic motion in this classical model in order to pick up all the classical information from a mixture of 
description. Our study is not to give differences between quantum dynamics and classical dynamics because different 
descriptions will definitely give different results such as using different phase-space distribution functions 24]. The 
only thing to do, in this paper, is to find out all the classical dynamics of atomic walking within a classical framework 
and make a try to show whether some nonclassical motions are more a classical collective behavior (or a long time 
statistical behavior) than a quantum wave description of an individual particle. In order to give a full picture of the 
classical motion for a quantum system easy for check, we closely investigate the classical walking of a single atom in 
a controlled optical lattice. We first derive the simple classical dynamic model in section UH and find out a variety 
of atomic walking behaviors of the model in section [TTTl Some approximate analytical solutions on specific conditions 
are given to understand the complicated behaviors of atomic walking in this unsolvable system. Section IIVI provides 
the simple discussions and the main conclusions. 
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II. THE SIMPLE CLASSICAL MODEL 



What we start with is a widely verified quantum model which describes the dipole interaction between an atom 
and a classical field named Rabi model |25| . The dynamic of a two- level atom with a dipole moment d interacting 
with a quasi-monochromatic plane-wave field, 



E(r,i) =ex#(r,t)cos(^-k-r), (1) 



is described by a Hamiltonian of 



p2 1 

H = -i— + -Mj q & z - Ml (r, t) cos (vt - k ■ r) (a~+(T + ) , (2) 
2m 2 

where the Rabi frequency is defined by fi(r,i) = ^ (r,t) ■ d/H. In Eq.([T]), S (r,t) is the temporal envelop of the 
wave field which can be definitely controlled in the experiment. The modulated Rabi frequency 0(r, t) expresses the 
coupling intensity of the field mode with the atomic dipole d, and & + is the atomic level raising (lowering) 

operator. Hamiltonian @ is a widely used model to study the interaction between a two-level atom and a classical 
field with different controlling modes denoted by 0(r, t). 

However, what we consider here is a simple case: a single atom interacts with a standing wave mode in a microcavity 
with only its axial direction, x, getting involved. In this case, the Hamiltonian becomes 

p 2 1 

H = 1 — huo(T z — Ml (x, t) cos vt (<T - +cr + ) , (3) 

2m 2 

where (x,t) depends on the axial cavity-mode with an envelop modulated by the input field. In a resonant mode 
field, Eq.® reduces to a so-called double resonance model [26| with an effective Hamiltonian of (see Appendix A for 
details) 

P 

H = hX cos(xi) cos (kx) , (4) 

2m 

where the key parameter x is the modulation frequenc y o f the field amplitude. Eq.Q is an extensively considered 
model [27h3Q| | along with the phase-modulated model (3l[ in quantum chaos. Here the Hamiltonian describes the 
dynamics of an atom (or a polar molecule) moving in an amplitude-modulated standing wave under an effective 
coupling constant A. However in quantum theory, this model is time-dependent and difficult to be solved in a normal 
technique [32J. For one single atom, we can introduce the following dimensionless classical variables H = (H)/Hu> r , 
p = (p)/hk, x — k{x), and dimensionless parameters A' = \/ui r , lj = x/^ri t' = ^rt, scaled by recoil frequency of a 
single photon uj r — hk 2 Jm. The label (O) = O means the expectation value of the corresponding quantum operator 
O and this is a very good approximation for one single atom. In atomic scale, a single atom can be treated as a 
classical point with its position identified by the center of mass, and no fluctuations induced by spatial superposition 
exist here. Therefore the classical model can be obtained (primes are omitted) by 

o 

P 

H(t) = — — Acosxcoswi. (5) 

According to Eq.([5]), the classical dynamic equations are 

9H dH , . 

x = —— = p, p— — - — — — Asinxcoswt, (6) 
op ox 



x + Xcos uit sin a; = 0. (7) 

Eq.Q looks simple, but it is a time-dependent nonlinear differential equation and no closed analytical solution for an 
arbitrary parameter is available (33| . 

The steady solutions of Eq.([7]) are clear at the antinode sites for x* = 0, ±7r, ±2ir, ■ ■ ■ and p* = 0, but they are not 
always stable because the trapping potential is time-dependent in the way of 

V (x, t) = — A cos {tut) cos a; = — ^ [cos (x + tot) + cos [x — tot)] , 
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FIG. 1: The effective potential along lattice x at different time t = 0.5 (solid line) and t = 1.8 (dashed line) with parameters 
A = l,u = 1. The time t is scaled by recoil frequency u) r and the position x is scaled by the wavelength of the cavity mode. 



which acts on the atom as a superposition of two potential waves propagating in opposite directions. In a high 
fineness cavity, this potential is usually used to control atomic motion and can be easily generated by a modulation of 
standing cavity-mode. FigQ] shows the varying potentials at two different times along the axial direction. The nodes 
(indicated by arrows) or antinodes of the potential lattice are fixed with a spatial period of 7r but the amplitude of 
the lattice wave changes with time with a period of t = 1~k /lu. The amplitude A at antinodes denotes the atom-lattice 
coupling indicating a trapping ability of the atom at antinode, so we can call it as a lattice depth. All the above 
lattice parameters can influence the dynamic of the atom and all can be controlled by an external optical field. 

However, due to atomic motion, the potential lattice felt by the atom also depends on the momentum of the atom. 
If we suppose that the velocity of the atom is a slowly varying quantity during a time period of [0,i], we can set 

x (t) = / p{r)dT = (p) • t w po ■ t, 
Jo 

where (p) is a time average of momentum which can be replaced by its initial value of pq. In this case, the optical 
force will be 

F(t) = -^- = - ^ [sin (w + po) t- Bfai (w-po)*], (8) 

which clearly indicates that there are two driven forces exerted on the atom with two different frequencies depending 
on atomic momentum, suggesting two nonlinear resonance points of the dynamics at to ± po [26j j . Therefore, the 
motion of the atom across the varying lattice field is dependent on two aspects: the lattice field and the state of the 
atom. Although this model described by Eq.([S]) is simple, it can display most of the dynamic behaviors found in 
[l9l - l23T | without directly including the internal dynamics. 



III. DIVERSITY OF WALKING BEHAVIOR 

A. Oscillation with random step-jumping 

Under the condition that the amplitude of the lattice is varying slowly, the cold atom (p < 2-\/A) will be trapped 
at the bottom of the lattice potential for a long time, oscillating around the antinodes at sinir* = 0, i.e. x* = 
±n7r, n = 0, 1, 2 • • • . In Fig(3J a typical walking process of this case is shown by a numerical simulation on Eq.©. 
Fig|2ja) displays a temporal position (black line) and momentum (gray line) of the atomic walking and Figl^b) is 
its orbit in the phase space. We can see that the atom first oscillates around one of antinodes (indicated by the 
horizonal dashed lines) and then jumps randomly in an integer steps of tt to a left or right antinode after a time 
interval of tt/lo. The random atomic jumping to a new site is due to the loss of stability of the former site when the 
coupling intensity X(t) = Xcos(ujt) turns from positive to negative. This classical jumping behavior demonstrated 
in this model will reversely affect the lattice field and can be traced by the transmission field from the cavity due 
to a motion-dependent detuning effect [U, [35[ . Although a similar dynamic behavior was revealed by Domokos and 
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FIG. 2: A typical walking trajectory of a cold atom in slowly varying standing wave with A = 1, uj — 0.02, the initial position 
and momentum xq — 0, po — 0.02. (a) The evolution of the atomic position (black line) and momentum (gray line), (b) The 
corresponding orbit in phase space. 



Ritsch [36(, the mechanisms of two systems are different. Our system is only for one single atom and no random 
Langevin-type noise is necessary during the dynamics. In order to understand more about this classical trapping and 
jumping behavior, we can analyze it by the following extreme approximations. 



1. Harmonic Oscillation with random jumping 



For the small modulation frequency of a; <C 1, the approximation of coswi ~ 1 is valid for a short time. The 
linearization for a cold atom (p <C vA for tightly trapping) around the stable positions of sins* = 0, such as sin a; ~ x 
around (for other site we can set x' = x — x*) can be used here. Then Eq.([7J will reduce to 

x + Xx^0, w<l, (9) 

where we set A > and this enables Eq.® a harmonic oscillation solution of 

i(t)»i (0) cos (Vxt) + ^ sin (VaY) , (10) 

which describes a main properties of atomic walking in this case. Eq. (|10[) indicates that the oscillation frequency in 
FigJ^a) is determined by the depth of the trap, A, and gives an important characteristic time of oscillation period, 
2tt/VX. As the temporal depth X(t) = X cos ujt changes slowly with time, the oscillation frequency will adiabatically 
follow with yj X(t). When X(t) becomes negative, the oscillation frequency yj X(t) — > iy/\X(t)\, and the atom will 
escape exponentially from the former site in a way of 

x(t)mx (0) cosh (Vxtj + sinh (yXtj . (11) 

Therefore the atom will conduct an oscillation followed by an escaping jump from the former site. In order to see 
the details of this motion, FigJ^a) exhibits a zoomed in orbit around one antinode in phase space. The black line 
is the real orbit and the thick gray ellipse stands for the harmonic oscillation of Eq. (lT0|) . where the potential lattices 
at time t = and at the end of one oscillation t — 27tt are depicted by the dashed lines for reference. As a decrease 
of X(t), the elliptic orbit in the phase space expands along x until X(t) decreases to a value that the atom can pass 
through it to another site at a location determined by the escaping momentum. Therefore the trapping time around 
antinodes for the cold atom is about tt/u as indicated in Figj2](a), which is estimated by the time when X(t) changes 
from positive value to negative value. In Fig|3Ja), the oscillating time of the atom around the starting site is about 
t = it/2lo = 257T, which is just the time taking by the trap depth A(i) decreasing from maximum to zero. 

Certainly, for a more rigorous approximation in above case of ut <C 1, the time dependent part will be cos(wt) f=a 
1 — (ujt) 2 /2, and the atomic motion can be better described by a parabolic cylinder function satisfying differential 
equation of 

x + X(l-^-)x = 0, (12) 
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FIG. 3: The enlarged orbit of a cold atom in phase space with uj — 0.02 and A = 1. The initial position and momentum are 
xo = 0.6, po = 0.1. (a) The orbit of the rigid simulation (solid black line) with harmonic approximation (gray thick line); (b) 
The approximation orbit of Airy solution of Eq, (|13|) . 



or even by an Airy function determined by 

x + - wt)x = 0, (13) 

where the time approximation cos(wi) = sin(7r/2 — cot) rj 7r/2 — cot is used. FigJ3Jb) demonstrates a better orbit of 
the Airy solution than a harmonic oscillation at a short time. But for a longer time the dynamic will dramatically 
deviate from the real orbit because of a divergence of Airy function. 



2. Pendulum rotation with random step-shifting 



When varying frequency to increases near to the oscillation frequency u> ~ vA/2, the trapping time tt/lo of the 
cold atom will decrease and the atom will quickly shift from one site to another like a pendulum swinging or rotating 
around different equilibrium sites just for a few periods. As the motion of atom covers a large range of x relative to 
stable antinodes in this case, the spatial linearized Eq.© will be invalid. Therefore Eq.© should be 

x + A sin x = 0, 

which is the well-known pendulum equation with A > 0. Actually, this situation corresponds to an atom moving in a 




FIG. 4: The phase space picture of a cold atom walking in a slowly varying field with A = 1, U) = 0.1 and the initial atomic 
momentum p(0) = 0. (a) The rigorous solution for initial position being x(0) — 0.4; (b) The pendulum solution of Eq. (|14[l by 
using a random step-shifting (step unit tt) position of Gaussian distribution with mean value of 0.4 and variance of n/2. 



stationary one-dimensional lattice. When A > 0, the integral solution of above equation takes the form of 

x 



x (t) — 2 arcsin 



siri Z1 S N fV\t + K (sin ^) , sin ^) 



(14) 
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with the initial conditions of x (0) = xq, x (0) = 0. The label SN(-, •) is the Jacobi elliptic function and K (•) is the 
complete elliptic integral of the first kind defined by the elliptic integral of 



K{k) 



*' 2 d9 



\j\-k 2 sin 2 i 



Walking orbits in the optical lattice with a higher modulation frequency lu = 0.1 are demonstrated in FigJU FigfJJa) 
is a rigid simulation and Figf^b) presents a pendulum samples of phase picture with random initial positions. We 
can see the long-time atomic walking can be somehow (not exactly) explained by the pendulum solutions with 
shifting equilibrium positions. As effective A(i) is actually time-dependent, the length of the pendulum will change 
gradually from positive to negative. When X(t) < 0, the stable site at x* — 0, ±2tt, zLAtt, ■ ■ ■ for X(t) > will become 
unstable and shift to a new stable site at x* — ±7r, ±37r, • • • , which is clearly shown in Figj^a) with a sequence of 
— 7r, 2ir, 37T, 2ir, 37T, Air , 7r, Oir , • • • . In FigfJJb) the thick black lines are pendulum solutions of Eq. (fT4| for A > and the 
thin dashed lines are for A < 0. We use a Gaussian random number to pick up initial positions for the pendulum 
solutions and the results reveal that a long time atomic walking in FigUJa) can be illustrated by Eq. (|14[) with a random 
shifting of equilibrium position of pendulum rotation. The simulation indicates that this behavior is a combination 
of pendulum oscillation with a random shifting of the optical lattice. 

According to above analysis, we can conclude that, in a slowly modulating optical lattice, the cold atom will feel a 
gradually changing wave trap, resulting in a harmonic or pendulum rotation followed by a random step-shifting. In 
this process, two characteristic times are important, trapping time 7r/w and oscillation period 2-k/\/X. In addition, 
the step-shifting between antinodes of the lattice has a random characteristic and can not be predicted because the 
jumping process is very sensitive to position and momentum. In this respect, the atomic walking in modulated 
standing wave is very similar to one-dimensional billiard problem that a minor deviation of position (or momentum) 
will leads to a dramatically different walking orbit. However, above conclusion is on the conditions that the optical 
lattice is slowly modulated and the atom is very cold. In the following parts we will find some other specific walking 
behaviors beyond these constraints. 



B. Random walking and Chaotic transportation 

When the modulation frequency of the lattice is further increased, the trapping time of the cold atom, tt/uj, will 
decrease to an extent that it is comparable to the shifting time between different sites, Tr/(p) ((p) is average momentum 
during a walking interval). In this case, the atom will exhibit a random like walking that is shown in FigJSJa) by 
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FIG. 5: (a) A typical random walking of a cold atom in varying standing field of u = 0.8, A = 2 with a;(0) = 0, p(0) = 0.4. The 
dashed line is the average motion of the atom, (b) The corresponding power spectrum of the random walking (solid) and the 
average motion (dashed). The inset picture is the pase-space distribution of the walking mapped to [0,27r]. 

a typical sample of moving sequence along optical lattice, displaying an unpredictable position of the atom in the 
periodically varying lattice. This random characteristic can be further verified by a power spectrum of this walking 
depicted in Figj5](b). The dashed line in FiglSJa) corresponds to the effective uniform motion of this walking with an 
average momentum of (p(t))x ~ 0.4 + 0.06 (T is the total walking time), whose fourier spectrum is depicted by a thick 
dashed line for reference. The inset picture of FigJSJb) is a long time pase-space distribution of the atom by mapping 
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its walking position into the first periodic region of the lattice [0, 2%]. The uniform phase-space distribution is a 
statistical proof of the random properties coming up in this deterministic system where regular and chaotic regions 
coexist in this walking process. 




0.0 0.5 1.0 1.5 2.0 1 2 3 4 

Coupling A Frequency co 

FIG. 6: (a) The Lyapunov exponent of the atomic walking versus coupling intensity A with uj = 0.8. (b) The Lyapunov 
exponent versus modulation frequency w with A = 2. The initial atomic position x(0) = 0, momentum p(0) = 0.4. 

In the sense of walking sensitivity to the atomic state, above behavior can be treated as a chaotic transportation in 
the standing wave lattice similar to the results in fl9l l23j . The estimation of Lyapunov exponents Xi can give a rough 
analysis of this chaotic property. This system is equivalent to a 3-dimensional autonomous system with one Lyapunov 
exponent being zero (the third dimension is 8 = ut). Because the sum of the exponents for a Hamiltonian system is 
null, i.e., J2i \ — 0> therefore the Lyapunov exponents of this system must be — r, 0, +r. Fig|5]gives a rough estimate 
of the largest Lyapunov exponent r changing with the trapping depth A (FigJBJa)) and with the modulation frequency 
u (FigjSJb)), indicating a weak chaotic transportation along the standing wave lattice. The chaotic transportation of 
this system under the influence of internal transition has been investigated by Argonov and Prants 0, [2(], [23| , but 
our results reveal that only the external walking of an atom in a slowly varying optical lattice can present the similar 
chaotic behavior. Certainly, the similar fractal tunneling time [37j through a lattice with a certain length can also be 
found in this walking. 



C. Quasi-periodic trapping state and its stability 

In some parametric regions the atom will be totally trapped by the quickly varying standing wave. FigJT^a) gives 
an example of this case when the frequency reaches to u = 2.8u r (left black bowknot orbit). In this case the atom 
takes a periodic or quasi-periodic oscillation around a stable site near initial position x(0). In the following, we will 
give some detail analysis about this trapped walking behavior. 

For an arbitrary varying frequency w, Eq.([7]) corresponds to a pendulum with changing length. The normal method 
of solving this nonlinear differential equation is the linearization approach. For the trapped atom that is very near to 
a stable antinode of sin a;* = 0, for example, a small x near zero is since ~ x, Eq. (|TJ) will be reduced to 

x + (A cos ut) x = 0, (15) 

which is a linear equation with a periodic coefficient and can be solved by Floquet's theorem (38j . The linearized 
equation (|15[) is a Mathieu equation 

d 2 x 

—j + [a - 2q cos (2t)) x = 

with a standard form of solution 

x(t) = e^Xit), 

where /i, in general, is a complex function of a and q called characteristic exponent, and X(t) is a periodic function. 
Explicitly, the solution of equation (TT51) reads 

x (t) = x(0)C (o, - g, f) + p(0)S (o, - g, ^ , (16) 
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FIG. 7: (a)The phase portrait of strict solution of Eq.(J7| (left black orbit) and the Mathieu solution of Eq. (|16p (right gray orbit) 
with parameters A = l,ui — 2.8 and the initial position and momentum being x(0) = 0.2, p(0) = 0.18. (b)The corresponding 
power spectrum of the atomic momentum for the strict solution (black line) and the Mathieu solution (gray line). 



where C (0, — ^f) and S (0, — ^k, ^f) are Mathieu even and odd functions. Above solution can be used to estimate 
a variety of behavior of Eq.([7| under the trapping condition (for that the linearization condition is valid). FigJTJa) is 
a comparison of Methieu solution (right gray orbit) with the strict simulation (left black orbit) in phase space when 
Methieu function is in its stable parametric region. The figures show a nice match of the linear approximation with 
Eq.© and the atom conducts a quasi-periodic oscillation around x* = 0, ±7r, ±27r, • • • . The quasi-periodic properties 
of this walking can be verified by the power spectrum of atomic momentum in FigJT^b), where several frequency 
components are manifest but the atom never exactly repeat its walking orbit in the phase space. Certainly, on some 
parameters, the periodic walking is also available. This can be discussed by the characteristic exponent of Mathieu 
function. The characteristic exponent, /j,(a, q) — /i(0, — 2A/oj 2 ), predicts that only for certain value of A and w can the 
solution be periodic. For parametric values of 2A/w 2 > 0.91, fi will be a complex number and Eq. fToT) is unstable (see 
the borderlines 2A/w 2 « ±0.91 shown in FigfSja)). Therefore we can determine the stable or unstable parametric 
region of the linearization equation under independent initial conditions. However in the real system, the modulating 
amplitude will include a small constant part £q [39j that is governed by 

d 2 x 

-jjj + (<?o + A cos ut) x = 0, 
and this enables a normal Mathieu solution of 

x (t) = x(0)C ^ j + P (0)S (-g-, --j, -j , (17) 

Fig J5] gives a stability diagram of solution Eq. fTT)) near stable positions x* with the shaded area indicating stable 
region of the trapping states. FigJSJ shows an irregular stable parametric region of the walking determined by the 
lattice depth A and the modulation frequency uj, with a symmetric structure for A and a fractal boundary for small u>. 
FigJSJb) also indicates that although the symmetric gap (indicated by the arrows in Fig(5Jb)) of the stable region, 2So, 
in small frequency u> region goes to zero when §q — >■ 0, the complicated structure of stable region will not disappear 
in FigJFJa) for a real walking. This result reveals that the dynamic stability of atomic walking around antinodes in 
a varying lattice is very sensitive to the fluctuation of the optical field [30]. Actually the fluctuation of the lattice 
field will definitely make the atomic walking complicated in the small modulation frequency region. Strangely, when 
modulation frequency is small, the atomic walking around antinodes can stay unstable no matter how large the depth 
of the lattice is. When the frequency becomes higher, the stable region will expand and the walking is more preferred 
to a stable motion. According to Floquet's theory, the characteristic exponent fj, of the trapping solution will be 
lim^-joo /i(0, —2X/lo 2 ) — 0, which indicates a free atomic walking limit in a quickly varying lattice. 

We should note here that the stable region of the real system Eq.([7]) is not only vulnerable to the control parameters 
A,w but also heavily depends on the dynamic itself, i.e., the atom's position and momentum are also key parameters 
to determine the stability of atomic walking [38] . The nonlinear position-dependent of walking in Eq. will introduce 
an effective part of £q to shift the atom along the lattice, which actually breaks the spatial symmetry of the lattice 
walking. Naturally, an increase of atomic momentum will definitely reduce the stability region of the real trapping 
solution because the heated atom will more easily escape from the trapping region. The stable region of the linear 
solution indicates that the atom can surely be trapped in a shallow lattice field (small A) with a higher modulation 
frequency u>. 
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FIG. 8: The stability region for trapping state of linear solutions Eq. (|17[) in the A — cj plane (a) for So — and (b) for So = 0.1 
The shaded area stands for the stable region. 



D. Free ballistic flying 



Finally, we consider two extreme conditions that the modulation frequency of the lattice is high {lu \/A) and the 
velocity of the atom is quick (p 3> 2vA). For a large modulating frequency uj, the effective coupling constant meets 
A cos(wi) — > and the equation can be simplified by 



0, u » 1, 



(18) 



which means a negligible influence of light on atomic walking and implies a free flying of the atom in optical lattice, 
i.e., x (t) x (0) + p(0)t and p(t) « p(0). In the other hand, when the speed of the atom is high, the atom also 
feels a quick oscillation of the force with a frequency of p + uj in one direction and p — uj in the opposite (see Eq.([8|)). 
Both above cases are displayed in Figj9](a) and Figj9jb) by direct simulations of Eq.([7]) where the thick black line, 
x(t), indicate a nearly free flying motion and its slope is roughly determined by the initial momentum po = p(0), 
showing a small relative change of atomic momentum during the walking process. The enlarged figures of the momenta 
between time interval [0, 30] shown by the insets of Fig|9] demonstrate a beat oscillation of the atomic momentum 
around the initial value, with the average momentum satisfying (p) po- Therefore, the oscillation behavior can be 
approximately described by Eq.@, 



dp 



— = — A sin (pot) cos (cot) . 



and an integral of it gives 



| A cos [(po - u)t] | A cos [(po + u)t] 

P\t) _ P0+ n + r, : 

2 po - oj 2 po + uj 



(19) 



which clearly reveals a beat oscillation of the momentum. Comparisons of solution Eq. (|19p (the gray lines of the insets) 
with the strict solution (the black lines of the insets) are depicted in Figj9ja) and Figj9jb). The simulation indicates 
a good agreement of Eq.(fT9"]) with the strict Eq.Q except for some minor differences. Besides, Eq. fHj)) also suggests 
a resonance behavior of the momentum when atomic momentum meets po ~ ±w. Under resonant condition, Eq. (|19[) 
becomes a bad approximation because the momentum will be divergent. However, as long as the condition of lj ^f~\ 
or po 2^f~\ is satisfied, the momentum enhancement by the resonant effect in real system remains small relative to 
the large initial momentum, and the atom still keeps its ballistic flying under resonance condition. However, if the 
initial atomic momentum is small, this resonance effect will be manifest only under small modulation frequencies (see 
Eq. ([T9"l) ) and this resonant behavior will also be suppressed by the nonlinear character of a real walking. Therefore, 
in conclusion, if the modulation frequency is high and the atomic momentum is large, the atom will almost conduct 
a free ballistic flying with its momentum taking a beat oscillation roughly around its initial value. 
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FIG. 9: The free flying of the atom (a) in a quickly varying field with ui = 10, p(0) = 0.1, a;(0) = 0.2, A = 1 and (b) in a high 
speed under parameters p(0) = 10, uj = 0.1, x(0) = 0.2, A = 1. The insets are atomic momenta zoomed in between a time 
interval of [0, 30] with the black lines for the strict equation and the grey ones for Eq. (|19[) , 



IV. CONCLUSIONS 



A single atom coupled with high-finesse cavities is a fundamental system for quantum research, and many theoretical 
and experimental works have been done on this system [l|-|5|, I34j - |36( ] . Directly based on quantum description, the 
quantum aspect of the dynamics are usually focused on, such as on nonclassical statistics, entanglement-induced 
effects, quantum information processing, which, contrarily, makes the classical contribution obscure. While some 
authors on this system concern more about the transmission of the light through a cavity to detect or control the 
atomic dynamics [34| . such as detecting atomic trapping states, or controlling atomic trajectories [13,|4l|. However 
all the atomic motion considered above are closely involved with the internal dynamic and the external dynamics is 
still mixing in atomic control. In order to investigate the external dynamic in a classical point of view, the internal 
variables should be decoupled or traced out. Raizen and his coworkers (42j | resort to a dynamical map (periodic kicked 
rotator) based on pulsed standing field to investigate the atomic dynamic behavior and find a good agreement with 
the classical dynamic in a noisy environment. Particularly, their experiment studies on the motion of cold cesium 
atoms in an amplitude-modulated standing wave of light |29l . l30l | have a very close relation to our study in this paper. 
However, we use a simple classical model under a resonant condition and find a more rich dynamic behaviors of one 
atom in the optical lattice beyond a direct influence of internal dynamics. 

As the fluctuating environment of the cavity mode, trapping a single atom in the cavity for a long time is turned 
out to be difficult. Therefore, sensitivity of the atomic motion to field modulation is a key problem for a practical 
one-atom control [3(| • In this paper, we closely consider the classical walking of an atom in an field-controlled standing 
wave and reveals a diverse dynamic region of atomic motion. In the parametric region of lower modulation frequency, 
an oscillation with classical random jumping is found for cold atom. With the increasing of frequency, random atomic 
motion, chaotic transportation and the quasi-periodic trapping state appear. If the modulation frequency is high 
or the velocity of the atom in the lattice is large, the atom will exhibits a ballistic fly with a momentum beating 
oscillation. The study of dynamic stability shows that the transition between these dynamical regimes is irregularly 
determined not only by the lattice field but also by the position and momentum of the atom. Our results indicates a 
wide parametric region of unstable walking and a susceptibility of stable motion to the field fluctuations of amplitude 
and frequency as well as to the motion itself. Therefore, this work gives a rich clue for atomic control in the cavity 
and also provides a useful insight into the dynamical behavior of atoms in a periodically varying lattice. 
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Appendix A: The derivation of Eq.(f4]) 

The dynamic of atom in a certain field can be manipulated by different configuration of E(r, t) according to the 
Hamiltonian of Eq.([2]). A standing wave, generated by an interference of two quasi- monochromatic field with an 
envelop, #(r, t), and the polarization, e^, propagating in opposite direction, will be 

E(r,i) = ie A <^ (r, t) [cos (vt + k • r) + cos (vt - k • r)] 
= B\§ (r, t) cos (k • r) cos (vt) , 

where v is the carrier frequency and the slowly varying temporal envelop, S (r, t), can be controlled in the experiment 
by changing the amplitude of the input field. If we stabilize the phase of the field in a linear cavity along x axis, the 
Hamiltonian of Eq.Q will be got. More generally, we can write down the Schrodinger equation in the standing wave 
as 



A 2 

^— + H a - Q (r, t) cos (vt) 
Zm 



^(r,t), 



where 

(r, t) = d • e\£ (r, t) cos (k • r) , 

and H a denotes the energy of the internal state: on ground state |1) with energy e\, on excited state |2) with energy 
£2- In real space, the atomic wavefunction can be written by a two-component form of 

* (r, t) = tp x (r, t) e- teit ' h \l) + <p 2 (r, t) e" 4£2 */ ri |2) . 

Under the rotating wave approximation, the Schrodinger equation for the two-component wave function is 

L(<PA-i( V R (r,t)e-^\(vi 
dt\f2j \n R (r,t)e iSt ^V 2 J 

where the detuning frequency S — ujq — v and the transition frequency ujq — (e 2 — ei) /K. The Rabi frequency is 
defined by 

V R (r, t) = = ^# (r, t) cos (k • r) . 

In the experiment, the carrier frequency v is adjusted to the transition frequency ujq (resonance 5 = 0), then 

' ' ' 1 2777 



dt \<P2 ) \-nn R {r,t) -|-v 2 ) \<f2 



In this case above equation can be decoupled by introducing [43 



ip+{r,t) = -i= [cpi (r, t) + <p 2 (r, t)] , V- (r, t) = [ip x (r,t) - ip 2 (r, t)} , 



with equations of 



d h 2 d h 2 

ih— ij+ = -—V 2 t() + -Ml R (r,t)^ + , ih—^ = -—V 2 ij_ + hn R (r,t)^-. 

Above equation for wave functions tp± (r, t) indicates that the atom on any state will feel two different potentials 

V± (r, t) = TtMn (r, t) = T—^S 1 (r, t) cos (k • r) . 

If the envelop function of the field is controlled by (amplitude modulation field) 

S (r, t) = g cos (*t) , 
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the effective potential will be 

d • G\ 

V (r,t) — =F5o — 2 — cos (xO cos (k • r) = S-A cos (xO cos (k • r) • 
For a linear cavity, the effective potential we pick up to consider is 

V (x, t) = h\ cos (xt) cos (kx) , 
where A is the effective coupling parameter and \ is the modulation frequency of the field amplitude. 
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